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We  have  measured  the  dissipation  and  frequency  of  nanocrystalline-diamond  nanomechanical 
resonators  with  resonant  frequencies  between  13.7  MHz  and  157.3  MHz,  over  a  temperature  range 
of  1.4-274  K.  Using  both  magnetomotive  network  analysis  and  a  time-domain  ring-down 
technique,  we  have  found  the  dissipation  in  this  material  to  have  a  temperature  dependence  roughly 
following  T^'^,  with  2^^«*10^"^  at  low  temperatures.  The  frequency  dependence  of  a  large 
dissipation  feature  at  —35-55  K  is  consistent  with  thermal  activation  over  a  0.02  eV  barrier  with  an 
attempt  frequency  of  10  GHz.  ©  2004  American  Institute  of  Physics.  [DOT  10.1063/1.1646213] 


Diamond  offers  the  largest  ratio  of  fETp,  among  mate¬ 
rials  used  for  nanomechanics,  where  E  is  Young’s  modulus 
and  p  is  density,  resulting  in  the  highest  resonant  frequency 
/o  for  a  given  device  geometry.  This  property,  combined  with 
high  thermal  conductivity  and  chemical  resistance,  makes 
diamond  an  attractive  material  for  nanoelectromechanical 
systems  (NEMS).  Recently,  Sekaric  et  aO  have  demon¬ 
strated  nanomechanical  resonators  in  nanocrystalline  dia¬ 
mond  at  room  temperature,  with  frequencies  of  6-640  MHz 
and  quality  factors  Q  of  2400-3500. 

In  addition  to  resonant  frequency,  dissipation  2^'  is  a 
critical  factor  in  NEMS  performance.  The  ultimate  perfor¬ 
mance  in  applications  such  as  microscopic  force  detection,^ 
signal  processing,^  and  fundamental  quantum  mechanical 
experiments,^'^  is  directly  related  to  the  dissipation.  To  illu¬ 
minate  the  possible  dissipation  mechanisms,  we  have  studied 
the  temperature  dependence  of  2  ^ '  for  four  diamond  reso¬ 
nators  with  frequencies  of  13.7  MHz,  55.1  MHz,  110.1  MHz, 
and  157.3  MHz,  between  1.4  K  and  274  K. 

Doubly  clamped  beams  were  fabricated  as  described  in 
Ref  1,  with  a  width  w  of  180  nm,  a  thickness  of  roughly  40 
nm,  and  lengths  f  =  2.5,  3,  4,  and  8  p.m.  A  35  nm  layer  of 
gold  was  deposited  on  top,  to  allow  magnetomotive  actuation 
and  detection®  in  fields  up  to  B  =  9  T.  The  sheet  resistance  of 
this  film  was  1.6  H/D  between  1.4  K  and  —20  K,  rising 
linearly  above  20  K  to  a  value  of  2.7  H/D  at  250  K.  The 
samples  were  placed  in  high  vacuum  (<  10^®  Torr)  on  a  He 
cryostat  with  a  base  temperature  of  1.4  K.  The  temperature 
was  monitored  using  a  carbon-glass  thermometer  and  regu¬ 
lated  using  a  resistive  heater,  both  mounted  on  the  copper 
sample  stage.  Given  the  large  temperature  dependence  of  the 
resonant  frequency  (Eig.  3),  temperature  stability  of  1  mK 
was  necessary  for  these  measurements. 

We  measured  the  quality  factor  and  resonant  frequency 
using  two  methods  that  rely  on  magnetomotive  effects.  Eirst, 
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we  used  conventional  network  analyzer  techniques.  Although 
this  approach  is  very  convenient  for  identifying  nanome¬ 
chanical  resonances,  extraction  of  accurate  quality  factors 
relies  on  a  detailed  model  of  the  embedding  and  stray  im¬ 
pedances.  In  order  to  verify  our  fitting  procedures,  we  have 
developed  an  alternate  time-domain  technique  which  is  not 
sensitive  to  these  impedances.  We  apply  a  10000  cycle  pulse 
train  at  the  resonant  frequency,  which  excites  the  mechanical 
oscillation  of  the  beam.  The  voltage  developed  across  the 
beam  is  amplified,  sampled,  and  stored  by  a  high-speed  stor¬ 
age  scope,  which  is  triggered  at  the  end  of  the  pulse  train. 
The  left-hand  inset  of  Eig.  1  shows  the  captured  free  decay, 
which  is  then  processed  using  a  fast  Eourier  transform  (EET). 
Eigure  1  shows  a  comparison  of  the  EET  power  spectrum  of 
the  free  decay  data  with  network  analysis  data.  Agreement 
between  the  two  techniques  for  the  extracted  quality  factor  is 
better  than  5%,  giving  us  confidence  in  our  network  analysis 
fits. 

Eigure  2  shows  the  central  result  of  this  letter:  The  dis- 


FIG.  1.  The  solid  line  shows  the  resonance  of  the  13.7  MHz  beam  at  T 
=  4  K,  measured  using  network  analysis.  The  filled  circles  show  the  results 
of  a  FFT  of  the  ring-down  transient;  the  raw  time-domain  data,  averaged 
500  times,  are  shown  in  the  left-hand  side  inset.  The  right-hand  side  inset 
shows  the  measured  dissipation  (open  circles)  as  a  function  of  magnetic 
field,  along  with  the  calculated  (solid  line)  magnetomotive  dissipation. 
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FIG.  2.  Dissipation  as  a  function  of  temperature  for  the  four  resonators. 
Data  for  the  13.7  MHz  resonator,  which  was  taken  at  5  =  9  T,  has  been 
adjusted  to  remove  the  effect  of  magnetomotive  damping,  using  Eq.  (1). 
(The  other  resonators  were  measured  at  a  low  enough  field  that  magneto¬ 
motive  damping  was  negligible.)  The  arrows  indicate  the  position  of  a  De¬ 
bye  peak;  the  inset  shows  an  Arrhenius  plot  of  \n(f)  vs  l/T  for  this  peak. 

sipation  2^'  of  the  diamond  resonators  as  a  function  of 
temperature.  In  order  to  understand  these  data,  we  have  con¬ 
sidered  a  variety  of  possible  dissipation  mechanisms:  Mag¬ 
netomotive  damping,^  thermoelastic  dissipation,^’^  mechani¬ 
cal  dissipation  in  the  gold  electrode, clamping  losses,"  and 
finally,  dissipation  in  the  diamond  itself.  We  will  address 
each  mechanism  below. 

Magnetomotive  damping,  which  results  from  electronic 
dissipation  in  resistive  elements  due  to  induced  eddy  cur¬ 
rents,  is  a  substantial  effect  for  these^  and  other*^  devices. 
This  mechanism  is  additive  to  the  “intrinsic”  dissipation 
2o  is  a  function  of  the  electromechanical  resistance  of 

the  beam,  Rg^={B()^QQ/2TTfQm,  and  the  impedance  of  the 
external  measuring  circuit,  Z^.^, : 

=  (1) 

where  B  is  magnetic  field  and  m  is  resonator  mass.^  Here, 

includes  the  series  combination  of  the  resistance  of  the 
gold  electrode  on  the  resonator  (—20-60  fl),  the  impedance 
of  the  1.5  m  semirigid  cable,  and  the  input  impedance  of  the 
preamplifier  (50  (1  or  3  kU).  The  right-hand  side  inset  in  Fig. 

1  shows  both  the  measured  dissipation  versus  magnetic  field 
and  a  fit  using  the  above  parameters,  with  the  cable  imped¬ 
ance  approximated  by  a  series  inductance  of  400  nH  and  a 
parallel  capacitance  of  110  pF.  Comparable  data  shown  by 
Mohanty  et  al.}^  which  was  attributed  to  the  presence  of 
charged  defects  or  electronic  donor  impurities,  can  be  ad¬ 
equately  fit  with  similar  parameters.  For  each  resonator  stud¬ 
ied  here,  we  have  measured  2^'  versus  B  and  determined 
the  largest  field  such  that  the  magnetomotive  damping  con¬ 
tributes  less  than  5%  of  the  total  dissipation. 

We  have  estimated  the  contribution  of  thermoelastic  loss 
and  find  this  to  be  negligible  for  all  of  the  resonators  pre¬ 
sented  here,  over  the  entire  temperature  range  studied.  This 
results  from  the  small  width  of  the  beam  (180  nm)  and  the 
high  thermal  conductivity  of  diamond,  which  combine  to 
make  the  thermal  time  constant  of  the  beam  several  orders  of 
magnitude  shorter  than  the  mechanical  period.* 

Another  contribution  to  the  dissipation  is  the  internal 
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losses  of  the  gold  film.  We  find  that  for  a  doubly  clamped 
beam  composed  of  two  materials,*"^ 

e^‘=Y^(emi+/3<2Au),  (2) 

where  /S=  fAuFAu/loia^Dia'^O.OS  in  our  devices.  Qj,  f,  ,  and 
are,  respectively,  the  quality  factor,  beam  thickness,  and 
Young’s  modulus  for  diamond  and  gold.  Using  Eq.  (2)  and 
the  measured  dissipation  in  thin  films  of  gold,*®  we  estimate 
that  the  contribution  of  the  gold  to  the  measured  dissipation 
is  less  than  5%  of  the  observed  dissipation. 

The  dissipation  of  the  shortest  (157.3  MHz,  € 
=  2.5  /rm)  beam  displays  an  additive  offset  compared  to  the 
other  three  beams.  This  offset  could  be  due  to  clamping  loss, 
which  is  expected  to  become  increasingly  important  in 
shorter  beams.  Calculations  of  energy  lost  to  semi-infinite 
two-dimensional  supports  for  cantilevers**  and  doubly 
clamped  beams*^  predict  dissipation  2ciamp“ 
where  /<=“  1 1.2  for  the  fundamental  in-plane  flexural  mode  of 
vibration.  For  the  2.5  /rm  beam,  2ciamp”“4.2X  10^*,  over  30 
times  greater  than  for  the  8  jjm  beam.  However,  this  estimate 
is  clearly  too  large  compared  to  the  observed  data,  which 
indicates  that  the  simple  two-dimensional  model  used  is  not 
sufficient  to  capture  the  behavior  of  the  beams.  Given  that 
this  loss  mechanism  is  predicted  to  be  increasingly  signifi¬ 
cant  for  high-frequency  resonators,  a  more  complete  model 
would  be  very  useful. 

Except  for  the  additive  offset  noted  above  in  the  157.3 
MHz  resonator,  the  dissipation  shown  in  Eig.  2  is  believed  to 
be  from  the  diamond  itself.  The  temperature  dependence, 
although  complex,  is  similar  to  that  observed  in  GaAs,** 
Si,*’*®  Si02,*"*  and  metallic  films,  such  as  Au  and  Al.*®  This 
apparent  universality  in  both  magnitude  and  temperature  de¬ 
pendence  provides  further  evidence  that  the  dissipation  arises 
from  a  broad  spectrum  of  tunneling  states,  similar  to  those 
found  in  amorphous  solids.*®  The  overall  trend  in  the  low- 
temperature  data  roughly  follows  a  power  law  with  an  expo¬ 
nent  of  —0.2,  which  has  yet  to  be  explained. 

A  prominent  feature  seen  in  this  data,  as  with  other 
materials,*’**’*®  is  a  Debye  peak  at  —35-55  K.  The  inset  of 
Eig.  2  shows  an  Arrhenius  plot  of  the  position  of  this  peak, 
which  displays  a  linear  dependence.  This  indicates  a  ther¬ 
mally  activated  process  with  energy  A£’«“0.02eV  and  an 
attempt  frequency  Fo^^lOGHz.  A  more  dramatic  peak  is 
observed  at  T=220  K  in  the  13.7  MHz  resonator.  It  is  diffi¬ 
cult  to  speculate  on  the  nature  of  this  peak,  as  it  was  not 
observed  in  the  other  beams  over  the  temperature  range  ac¬ 
cessible  by  our  apparatus. 

Eigure  3  shows  the  temperature  dependence  of  resonator 
frequency,  which  is  similar  in  magnitude  and  sign  to  shifts 
observed  in  Si  and  GaAs.**’*®  However,  the  increasing  fre¬ 
quency  at  higher  temperatures  (>160  K)  in  the  13.7  MHz 
and  the  110.1  MHz  resonators  is  different  behavior.  (The 
55.1  MHz  resonator  was  destroyed  before  we  were  able  to 
perform  measurements  at  higher  temperatures.)  Although  we 
suspect  this  effect  is  a  result  of  relative  thermal  expansion, 
calculating  the  size  of  thermal  strain  is  complicated  by  the 
complex  structure  of  the  samples:  Three  thin  films  (Au,  dia¬ 
mond,  and  Si02)  with  unknown  stresses,  all  atop  a  Si  sub¬ 
strate.  A  simple  model  incorporating  only  the  length  change 
to  AlP  license  or  copyright,  see  nttp://apl.aip.org/apl/copyright.jsp 
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FIG.  3.  The  shift  in  resonant  frequency  5///o  as  a  function  of  temperature 
for  the  four  resonators,  plotted  on  a  logarithmic  scale.  The  solid  and  dotted 
lines  show  the  predictions  of  a  model  based  on  the  thermal  expansion  of  Si 
and  diamond. 

and  tension  produced  in  the  resonator  by  the  thermal  expan¬ 
sion  of  the  Si  substrate  and  diamond  film  predicts  a  mini¬ 
mum  in  <5///o  around  140-160  K,  with  the  appropriate  order 
of  magnitude.  However,  the  deviation  at  low  temperatures 
indicates  that  factors  not  yet  understood  also  play  a  signifi¬ 
cant  role.  Interestingly,  the  magnitude  of  the  frequency  shift 
does  not  depend  monotonically  on  resonator  length,  which 
may  be  a  result  of  mesoscopic  variations  of  the  mechanical 
properties  from  sample  to  sample.  This  might  be  expected, 
since  the  diamond  grain  size  (—5-40  nm)  can  be  a  signifi¬ 
cant  fraction  of  the  resonator  width  (180  nm). 

In  summary,  we  have  measured  the  dissipation  of 
nanocrystalline-diamond  nanomechanical  resonators  be¬ 
tween  13.7  MHz  and  157.3  MHz.  We  find  that  the  magnitude 
and  temperature  dependence  of  2  ^ '  are  very  similar  to  those 
observed  in  other  materials,  suggesting  a  similar  mechanism 
for  dissipation.  The  power  law  — ^  observed  at  low 
temperatures  here  and  with  other  materials  begs  further  study 
and  understanding.  To  help  answer  these  questions,  measure- 
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ments  of  this  and  other  materials,  including  single-crystal 
diamond,  at  ultralow  temperatures  (10  mK-1  K)  would  be 
very  illuminating.  In  addition,  by  combining  the  time- 
domain  technique  demonstrated  here  with  very-high- 
sensitivity  readout  amplifiers,'^  it  may  also  be  possible  to 
observe  discrete  dissipation  events.  This  would  offer  further 
insight  into  the  mesoscopic  behavior  of  these  nanomechani¬ 
cal  devices. 
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